Astronomy & Astrophysics manuscript no. pap'rev's 


February 5, 2008 


(DOI: will be inserted by hand later) 





Stokes diagnostics of simulated solar magneto-convection 

S. Shelyag 1,2 , M. Schiissler 2 , S.K. Solanki 2 , and A. Vogler 2 



o 
o 



OS 



> 
o 

On 

m 
O 

o 

I 

o ■ 
a: 



1 Solar Physics and Space Plasma Research Centre, Department of Applied Mathematics, University of Sheffield, 
Hicks Building, Hounsfield Rd, S3 7RH, UK 

2 Max-Planck-Institut fur Sonnensystemforschung, 37f9f Katlenburg-Lindau, Germany 
February 5, 2008 

Abstract. We present results of synthetic spectro-polarimetric diagnostics of radiative MHD simulations of solar 
surface convection with magnetic fields. Stokes profiles of Zeeman-sensitive lines of neutral iron in the visible 
and infrared spectral ranges emerging from the simulated atmosphere have been calculated in order to study 
their relation to the relevant physical quantities and compare with observational results. We have analyzed the 
dependence of the Stokes-/ line strength and width as well as of the Stokes- F signal and asymmetries on the 
magnetic field strength. Furthermore, we have evaluated the correspondence between the actual velocities in 
the simulation with values determined from the Stokes-/ (Doppler shift of the centre of gravity) and Stokes-F 
profiles (zero-crossing shift). We confirm that the line weakening in strong magnetic fields results from a higher 
temperature (at equal optical depth) in the magnetic flux concentrations. We also confirm that considerable 
Stokes-U asymmetries originate in the peripheral parts of strong magnetic flux concentrations, where the line 
of sight cuts through the magnetopause of the expanding flux concentration into the surrounding convective 
donwflow. 
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1. Introduction 

Measurements of the photospheric magnetic fields of the 
Sun are mostly based on polar imetric observations of spec- 
tral l ines (for an overview, see lSolan ki 1993; Solanki et al 
200G). In particular, two line pairs of neutral iron (Fei), 
one in the visible wavelength range at 630.15 nm (effec- 
tive Lande factor g e ff — 1.66) and 630.25 nm (g = 2.5), 
and one in the infrared at 1564.85 nm (g = 3) and 
1565.28 nm (g e ff — 1-53) have often been used for po- 
larimetric studies since they are not blended by other 
lines, have large Lande factors, and because the small 
wavelength separation of each pair permi ts simultane- 
ous ob s ervatio n s with spectrographs (e.g., IStenflo et al 
1987bt iRabinl Il992l: ISolanki et alj Il992t i Lites et al 



1996t iMartmez Pillet et all Il997t ISigwarth et all Il999t 
Stolpe & Kneerl bOOfj ISocas-Navarro fc Litesl l200l . For 
magnetic field diagnostics in the visible spectral range, 
the Fei lines at 630.15 nm and 630.25 nm have the 
advantage of being less temperature sensitive than the 
other widely used l ine pai r 524 .7 nm and 525.0 nm 
(e.g.. IStenflol Il973l: Iwiehn Il978l : IStenflo et alj Il987a : 
Grossmann-Doerth et al.l ll996). Compared to the line pair 
in the visible, the infrared lines considered here are formed 
deeper in the atmosphere. Infrared spectral lines with suf- 
ficiently large Lande factors are well suited for the diag- 



nostics of weak magnetic fields since the Zeeman splitting 
grows proportionally to the square of the wavelength while 
the Doppler width increases only linearly. 

In this paper we carry spectro-polarimetric diagnos- 
tics of 3D rad iative MHD simulations obtained wit h the 
MURaM code (|Vogler et al.l2003H Vogler et al.l2005f) . con- 
centrating on the lines Fei 630.25 nm in the visible and 
Fei 1564.85 nm in the infrared. While the analysis of 
Khomenko et al. (|2005allrJ l considered mixed-polarity re- 
gions with zero net vertical flux through the computa- 
tional box, our study is based upon two simulations with 
non- vanishing net flux corresponding to horizontally aver- 
aged magnetic field strengths of 200 G and 10 G, respec- 
tively. These simulations are considered to represent small 
parts of active regions in early (plage) and late stages, re- 
spectively, of their development. We calculate synthetic 
profiles of the Stokes parameters for the visible and in- 
frared Fe I lines based upon the physical quantities in the 
simulated solar photosphere, analyze their relation to the 
atmospheric structure, and compare them with observa- 
tional data. In particular, we statistically analyze the ori- 
gin of asymmetries of Stokes V profiles originating from 
different regions in the computational box and relate them 
to the velocity and magnetic field profiles along the corre- 
sponding lines of sight. 
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The paper is organized as follows. Sect. 2 gives a brief 
description of the simulations underlying our analysis. The 
line synthesis on the basis of the simulation data is ex- 
plained in Sect. 3. The results are presented and discussed 
in Sect. 4, and Sect. 5 summarizes our conclusions. 



2. MHD simulations 



The MURaM code (|V6glerl 120031 : IVogler et all 120051) 



integrates the system of MHD equations on a three- 
dimensional, equidistant cartesian grid. The code is par- 
allelized using a domain decomposition scheme. The up- 
per boundary of the computational domain is assumed to 
be closed. It is located in the upper photosphere near the 
temperature minimum, where the density is rather low, so 
that the influence of the closed boundary on the dynam- 
ics of the granulation is negligible. The lower boundary 
is located in the convectively unstable layers of the up- 
per convection zone. An open boundary is implemented 
in order to allow free motions of the fluid through the 
boundary. The code includes non-grey (opacity binning) 
radiative transport and an equation of state taking into 
account partial ionization of the 11 most abundant chem- 
ical elements in the solar photosphere. 

The size of the computational domain for the simula- 
tions considered here is 6000 x 6000 x 1400 km 3 with a 
grid size of 288 x 288 x 100 grid cells. The corresponding 
resolution is 20.8 km in the horizontal directions and 14 
km in the vertical. The simulation starts with a plane- 
parallel atmosphere extending between 800 km below and 
600 km above the average level of continuum optical depth 
unity at 500 nm. After convection has developed and both 
the outgoing energy flux and the total kinetic energy have 
reached stationary values (with only short-term fluctua- 
tions on the granulation time scale), a homogeneous ver- 
tical magnetic field of 200 G and 10 G, respectively, has 
been imposed. After the decay of all transients caused by 
the introduction of the magnetic field, we have continued 
the simulations for about another hour solar time to be 
sure that the results have become independent of the ini- 
tial distribution of magnetic flux. 

The snapshot from the 200 G run (representing a small 
part of a plage region) is taken about 160 minutes after 
the start of simulation and 95 minutes after the verti- 
cal magnetic field was introduced. The snapshot from the 
10 G run (late decay phase of an active region, or base 
of a coronal hole) is taken about 165 minutes after the 
start of the simulation and 45 minutes after the vertical 
magnetic field was introduced. Maps of the vertical com- 
ponents of magnetic field and flow velocity on the surface 
7500 = 1 f° r the two simulation snapshots are shown in 
Fig. [T] They illustrate the concentration of the magnetic 
flux in the intergranular downflow lanes and the forma- 
tions of flux sheets (mainly in the 200 G snapshot) and 
tube-like magnetic structures. 



3. Line synthesis 

We have calculated Stokes profiles for the Fei lines at 
630.15 nm, 630.25 nm, 1564.85 nm and 1565.29 nm. The 
main properties of the transitions corresponding to these 
lines, such as the central wavelength (Ao), the transi- 
tion term, the effective Lande factor (g e g), the excita- 
tion potentials of the lower level (x e )j and the weighted 
logarithmic oscillator strengths (logg*/) are given in 
Table [TJ We have us ed a value of e = 7.43 fo r the 
iron abundance (e.g.. [Shchukina fc Truiillo Bueno 



2001 



Bellot Rubio fc Borreroll2002t lAsplund et al.ll2000bl ). The 
log(q *,f) values for the in frared lines have been taken 
from iBorrero et al. ( 20031 ) . For the 6 3 0.15 nm line we 



use the 
630.25 r 



value given by iBard et al 
m we use the same 



value 



thors (e. g. ICab rcra Solana et alJ 120051 : iKhomenko et al 



2005al ibl: IBorrero et allT 2006) 



( 19911 ) while for 
as previous a u- 



Ao [nm] 


Transition 




Xo [eV] 


log(ff7) 


630.15 


z »p» - e°D 2 


1.66 


3.654 


-0.718 


630.25 


2 5p0 _ e B Do 


2.5 


3.686 


-1.235 


1564.85 


e 7 D 1 - 3d 6 4s5p 7 L>? 


3 


5.43 


-0.675 


1565.29 


f T D 5 ~ (9/2)[7/2]2 


1.53 


6.25 


-0.043 



Table 1. Parameters of the spectral lines considered in 
this study. 



The line profile calculations have been carri ed out by 
STOPRO routine (|Solankilll987t iFrutigerl l2000h . It com- 
putes the full Stokes vector by integrating the radiative 
transfer equations under the assumption of local thermo- 
dynamic equilibrium (LTE), including Zecman splitting 
and magneto-optical effects. In order to obtain quantities 
required for line profile calculations which are not deliv- 
ered directly by the MURaM result files, the MODCON 
routine was used. It determines the LTE ionization equi- 
librium for a given chemical composition and derives elec- 
tron pressure, continuum optical depth, and continuum 
absorption coefficient. 

We have developed the program LINE, which calcu- 
lates the Stokes parameters for the 3D atmosphere simu- 
lated by the MURaM code. The program uses STOPRO 
and MODCON as subroutines. The program was paral- 
lelized for computational efficiency of the line profile cal- 
culations in large 3D computational domains. 



4. Results of Stokes diagnostics 

4.1. General properties 

We have calculated Stokes profiles for both pairs of Fei 
lines, taking vertical lines of sight through each of the 
288 x 288 horizontal resolution elements of the computa- 
tional box. 

In order to demonstrate the general consistency of syn- 
thetic line profiles from simulations with observations, we 
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show in Figure [5] a comparison of the spatially averaged 
Stokes-/ profiles of the Fe I line pairs in the infrared (upper 
panel) and in the (lower panel) from the (B z ) — 10 G sim- 
ulation (diamond symbols) with the obs erved spectrum 



of the quiet Sun (jDelbouille et alj 1 19731 . black curves) 



Observed and simulated line profiles agree to within a 
few percent. A more detailed comparison of mean pro- 
files would require considering a large number of snapshots 
from the simulation run covering several periods of the 5- 
minute p-mode oscil lation, which is beyon d the scope of 
this paper (see, e.g.. lAsplund et al.ll2000al ). 

In what follows, we shall concentrate on the results for 
one line in the visible (630.25 nm, for simpler notation 
abbreviated to 630.2 nm in what follows) and one line in 
the infrared (1564.85 nm, abbreviated to 1564.8 nm). We 
define a number of useful quantities to describe the prop- 
erties of the synthetic profiles, namely, the line strength, 

1_ f(Ic-I) 
A J I c 

where I is the (wavelength-dependent) intensity in the 
line, I c the continuum intensity, and Ao the wavelength 
of the line center in the rest frame. It is clear from the 
definition that the line strength is equal to the equivalent 
width normalized by the central wavelength. Furthermore, 
we consider the unsigned Stokes-V area, defined as 



S= — 



-dA, 



(1) 



A v = 



— dA, 



(2) 



where V is the Stokes V parameter, and the Stokes-V area 
asymmetr^ 



5A V = 



JVdX 



I\V\dX' 

as well as the Stokes-V amplitude asymmetry, defined as 



(3) 



Sa v = flail - |a,.|)/(|a b | + |a r |), 



(4) 



where ah is the amplitude of the 'blue' lobe of the Stokes- 
V profile and a r is the corresponding quantity of the 'red' 
lobe. In all cases above, the integrals are taken over the 
full width of the line. 

Figure [3] shows some results obtained for the snap- 
shot from the (B z ) = 200 G (plage) run. The normalized 
continuum image for A = 630.2 nm shows localized bright- 
ness enhancements corresponding to the concentrations of 
strong vertical magnetic field. The latter are clearly de- 
tectable in the maps of the unsigned Stokes-V area (mid- 
dle left panel for the 630.2 nm line, bottom left panel 
for 1564.8 nm) as can be seen by comparing with the 
magnetic field map in the top left panel of Fig. [1] While 
the unsigned Stokes-V area in both lines follows the to- 
tal field strength rather well, the larger Zeeman sensitiv- 
ity of the 1564.8 nm infrared line is clearly reflected in 



1 The sign of the net magnetic flux through the simulation 
box has been chosen such that, for a normal two-lobed V- 
profile, 8 Ay is equal to the result of the usual definition — 
|A r |)/(|^46| + |-Ar|), where Ab and A r are the areas of the 'b lue' 
and 'red' lobes, respectively (e.g.. ISolanki fc Stenfldll984l h 



the maps, more than compensating for its smaller line 
strength. The larger Ay value of this line may partly also 
result from its lower formation height and from the fact 
that strong flux concentrations show a larger height gra- 
dient of the field strength (compare the magnetic field 
maps in Fig. [1] with those in Figs. [7] and |5J. The corre- 
sponding maps of the line strength (top middle and top 
right panels, respectively) reveal significant line weakening 
in the regions of strong magnetic field. The maps of the 
Stokes-V area asymmetry (central panel and bottom mid- 
dle panel, respectively) and of the amplitude asymmetry 
(middle right and bottom right panel, respectively) show 
strong asymmetries in the periphery of magnetic flux con- 
centrations, while their central parts exhibit almost no 
asymmetry. The predominantly positive values indicate 
that in most cases the 'blue' lobe the of Stokes-V profile 
has larger area and amplitude than the 'red' lobe. Lines 
of sight through the peripheral parts of a flux concentra- 
tion that expands with height cross its boundary and, in 
the presence of strong external flows, thus sample signifi- 
cant magnetic field and velocity gradien ts. As a result, the 



Stokes-V profiles become asymmetric (jllling et al.l 11975) 



with a stronger bl ue lobe in the case of dominati ng ex- 
ternal downflows ( Grossmann-Doerth et al. 1988b! ) . The 
asymmetry arises from the fact that the wavelength shifts 
due to the Doppler effect and the Zeeman splitting have 
the same sign for one lobe, but opposite sign for the other 
lobe. The effect is strongest when both shifts are of sim- 
ilar magnitude; in addition, the area asymmetry depends 
on lin e saturation ( Grossmann-Doerth et alJll989t Solanki 
1989). This explains why the asymmetries are consider- 
ably larger for the visible line, which is stronger (more sat- 
urated) and has a smaller Zeeman splitting, better match- 
ing the Doppler shift. On the other hand, the spatial distri- 
butions of 5 Ay and Say from both lines are very similar. 

Figure [4] shows the corresponding maps for the (B z ) = 
10 G run, which could represent an almost decayed unipo- 
lar part of an active region or the field underlying a coro- 
nal hole. The maps of unsigned Stokes-V area (middle 
left and bottom left panels) show a significant signal only 
in a few isolated patches of concentrated magnetic flux. 
The infrared line displays an enhanced Ay not only in 
the weak-field regions, but also in the stromg flux con- 
centrations. This does not imply that the Stokes-V am- 
plitude is also larger for this line, since the V lobes can 
become very broa d owing to the vertical field gradient 



(jZaver et al. 1989). The maps of the Stokes-V area and 



amplitude asymmetries are very complex and show almost 
no correspondence to the structures in the continuum im- 
age. In this case, the strongest asymmetries arise from very 
weak magnetic fields in granules, with some preference for 
the granule edges. No sign of the asymmetry parameters is 
preferred, with the 630.2 nm showing much larger values. 

Synthetic 'slit spectra' of all four Stokes parameters 
along the vertical yellow line on the continuum image for 
the 200 G snapshot (upper left panel of Fig. [3]), are shown 
in Figs. [5] and [6] The width of the artificial slit is equal 
to the horizontal cell size of the simulation (about 20 km 
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on the Sun). The spectra represent 'ideal' cases; for a di- 
rect comparison with actual observations, one would have 
to take into account the instrumental profile and finite 
spatial resolution. The synthetic spectral lines show the 
characteristic 'line wiggles' from the Doppler shifts due 
to the granular up- and downflows. The Stokes-1^ spectra 
(upper right panel) show strong signals where the slit cuts 
through an extended flux concentration (in the upper half) 
and through a narrow flux sheet (at ~ 1800 km). The I- 
profiles of the lines are weakened and split (more strongly 
so in the case of the infrared line) at these locations, while 
strong line wiggles in both Stokes-/ and Stokes- V indicate 
the presence of significant internal structure of the large 
magnetic flux concentration. Note that some of the varia- 
tion may be due to the fact that the large, elongated flux 
sheet is bent, so that central and peripheral regions are 
sampled at different locations along the slit. The spectra 
of Stokes-Q and U show rather weak signals, consistent 
with the lack of significant horizontal magnetic field com- 
ponents in the strong flux concentrations. As expected 
from its larger Zeeman sensitivity, the infrared line dis- 
plays larger values of Q and U. . 

In order to compare spectroscopically determined ve- 
locities with their actual values in the simulation, we de- 
fine the line-of-sight velocity, Av, according to the Doppler 
shift of the Stokes-/ profile, viz. 



much better representation of the height to which the 



At; 



c(A - A cg ) 
An 



(5) 



where c is the velocity of light, An is the wavelength of the 
line center in the rest frame, and A cg is the wavelength 
corresponding to the center of gravity of the line. The 
latter is defined as the first moment of the profile, 



Ar 



J(I C - I)\d\ 
Rlc-I)d\-' 



(6) 



where I c is the continuum intensity and I is the line in- 
tensity. 

Fig. □ (for the (B z ) = 200 G simulation) and Fig. 
(for (B z ) = 10 G) give a comparison between the velocity 
Av determined from the Doppler shift of the 630.2 nm 
(middle left panel) and 1564.8 nm lines (middle right 
panel), respectively, with the actual flow pattern of the 
simulation at the surface t^qo = 0.1 (top left panel). 
This level can be seen as a very rough representation 
of the optical depth range contributing to the spectral 
line shift. The maps show that the actual velocities 
at that height and the spectroscopically determined 
velocities qualitatively agree. The somewhat larger flow 
velocities returned by the infrared line are consistent 
with its lower height of formation since the flow speeds 
tend to increase with depth. However, the corresponding 
scatter plots given in Fig. [5] show that the scatter 
of the individual points is considerable. This is not 
surprising in view of the extended height range of line 
formation (or better, the finite width of the velocity 
response function in combination with significant vertical 
velocity gradients). In fact, it has been shown that a 



termined via resDonse functions (Beckers & Milkcv 


19751 


Landi Deal'Innoccnti & Landi Dee 


'Iinnocenti 


19771 


Grossmann-Doerth et al.l Il988a: 


Del Toro Iniesta et al. 


1994 ISanchez Almeida et al.l 1996; Cabrera Solana et al. 


20051: iDelToro Iniestal 20031 Sect. 


10.4.3). Nevertheless, 



it is interesting to see that, for a first guess, the level 
T 500 = 0.1 is not a completely unreasonable choice. 

The bottom panels of Figs. [7] and [8] show maps of the 
line width, AA, which is defined here as the standard de- 
viation of the profile (square root of the second moment) : 



AA = 



l J(I c -I)(X-X c ydX 

j(i c - i)dx ■ 



(7) 



The maps of AA (bottom left panel for 630.2 nm, bottom 
right panel for 1564.8 nm) reflect the different magnetic 
sensitivity of the visible and infrared Fe 1 lines. In the case 
of 630.2 nm, the line width is actually decreased in the 
bulk of the magnetic flux concentrations, showing that its 
Zeeman sensitivity is insufficient for a significant broaden- 
ing. On the other hand, there is a clear increase of AA in 
the periphery of the flux concentrations, where the field is 
weaker but strong downflows prevail. In connection with 
the fact that the line width is found to be also enhanced 
in non-magnetic intergranular downflow lanes (see Fig. [5]) , 
this indicates that the strong velocity gradients in down- 
flows are responsible for the broadenin g of the visible line 
( Nesis et al. Il996 ; Solanki et al.1 Il996l ) . In the case of the 



more Zeeman-sensitive 1564.8 nm infrared line, the maps 
show a clear correspondence between the magnetic field 
strength and the line width. Doppler broadening is rele- 
vant for this line as well, since some enhancement of the 
line width can be seen in all downflow lanes, even in the 
absence of a strong field. 

4.2. Cut through a magnetic feature 

A two-dimensional vertical cut through a sheet-like mag- 
netic feature along the inclined white line in the upper left 
panel of Fig. [3] is shown in Fig. [TU1 The four panels dis- 
play the distributions of temperature (top left), gas pres- 
sure (top right), magnetic field strength (bottom left), and 
vertical velocity (bottom right). The green line on the im- 
ages indicates the T500 = 1 level. The analyzed flux sheet is 
rather shallow, the flux concentration extending only to a 
depth of about 500 km below the average level of T500 = 1 ■ 
With respect to its surroundings at the same geometrical 
depth, the gas in the interior of the magnetic flux sheet 
is cooler and has a deficit in gas pressure, which is al- 
most completely compensated by the magnetic pressure, 
i n accordance with the approximation of thin flux tubes 
(IShelvadl2004h . There are strong downflows with velocities 
up to 3 kms -1 at the periphery of the magnetic flux sheet, 
while the flow speeds in its interior are smaller, exhibiting 
an upflow of up to 600 ms _1 near t^qq = 1. Note that 
the velocity plot in Fig. [T0l shows only the vertical compo- 
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nent: since the flow is largely subsonic, a divergence of the 
vertical mass flux at its stagnation points is balanced by 
a convergence of the horizontal mass flux, and vice versa. 
The height level of t 500 = 1 is depressed by about 150 km 
in the region of strong magnetic field. This is caused by 
the lower gas temperature and density in the flux sheet, 
an effect analogous to the Wilson depression in sunspots. 

Fig. [11] shows the profiles of magnetic field and veloc- 
ity along several vertical lines of sight (indicated by the 
short vertical lines at the top of the panels in Fig. [TO"]) 
together with the emergent Stokes-/ and Stokes-F pro- 
files at these locations. Each line of sight corresponds to 
one row of panels in Fig. [Til which show (from left to 
right) the line-of-sight components of the magnetic field 
and the velocity field as a function of t 5Q0 , and the Stokes- 
/ and V profiles of the 1564.8 nm and 630.2 nm lines. 
The first row in Fig. [TTJ represents the center of the flux 
sheet, showing a slow upflow and large field strength in 
the height range relevant for line formation, so that the 
Stokes profiles are split and nearly symmetric. The field 
strength drops with height in accordance with the horizon- 
tal pressure balance. Moving towards the periphery of the 
flux concentration, the lines of sight cut through the mag- 
netopause of the expanding magnetic structure into the 
surrounding downflow (cf. the last three rows of Fig. [TTJ 
which show an increase of the field strength with height in 
the lower photosphere). As a result, the Zeeman splitting 
is smaller and the ^-profiles, particularly for the visible 
lin e, become st ro ngly asymmetric . 



Tiling et al 



(J1975J) have shown that asymmetric 
Stokes-V profiles arise in the presence of magnetic field 
and velocity gradients along the line of sight, the sign of 
the asymmetry being determined by the sign of the prod- 
uct of both gradients. On this basis, a mechanism leading 
to strongly asymmetric (and unshifted) Stokes-V^ profiles 
in the periphery of a s tatic magnetic flux conce n tration 
has been su g gested by iGrossmann-Doerth et al ] (|l988ch 
and Solanki "(1989). These authors assumed an idealized 
two-component atmosphere model: the upper part of the 
atmosphere has a magnetic field and is static while the 
lower part is non-magnetic and exhibits a systematic ver- 
tical flow. The Stokes- profile formed along a vertical line 
of sight that cuts through both components then becomes 
asymmetric with an u nshifted zero-crossing wavelength. 
Solanki fc Pahlkd (|l988h showed that the blue wing of the 
Stokes-V^ profile is stronger than the red wing (positive 
area asymmetry) if the inequality 

dr dr 

holds, where B\ os and i>i os are the line-of-sight components 
of the magnetic field and the flow velocity, respectively. 
Note that positive values of the velocity correspond to up- 
flows. The asymmetry is positive in regions where the ve- 
locity and the magnetic field gradients have the same sign, 
such as in the periphery of the flux sheet shown in Figs.fTDl 
and QTJ The area asymmetries of the ^-profiles corre- 
sponding to the five lines of sight considered in Fig. [TTJ 



are in agreement with this rule (compare Table [21 with 
the first two columns of Fig. [TTJ. Note that the Stokes- 
V profiles are formed higher in the atmosphere near the 
periphery of the flux tube, since the they obtain a larger 
contribution from layers with a stronger field. The fact 
that the 630.2 nm line generally shows larger ^-profile 
asymmetry than 1564.8 nm line is due to the combination 
of st ronger saturation and weaker splitting of the visible 
line ( Grossmann-Doerth et al. 19891 : Solanki 1989). 





630.2 nm 


1564.8 nm 


l.o.s. # 


Sav 


SAv 


Sa 


8A 


1 


-0.02 


0.03 


-0.11 


0.01 


2 


0.06 


-0.07 


0.16 


-0.03 


3 


0.31 


0.05 


-0.12 


0.04 


4 


0.51 


0.43 


-0.02 


0.33 


5 


0.72 


0.59 


0.13 


0.27 



Table 2. Stokes-V^ amplitude asymmetry (Say) and area 
asymmetry (SAy) for the 630.2 nm and 1564.8 nm lines 
along 5 lines of sight (l.o.s.) in Figs. [TUl and [HI The l.o.s. 
# 1 is near the center of the flux sheet while l.o.s. # 5 is 
in its far periphery. 



4.3. Statistical properties 

We now consider the complete set of Stokes profiles cal- 
culated for vertical lines of sight through each of the 
288 2 horizontal resolution elements of 20.8 x 20.8 km 2 size 
each, which are provided by a simulation snapshot. We 
study the dependence of their properties on the mag- 
netic field and flow structure using the snapshot from the 
(B z ) = 200 G simulation. 

Scatter plots of the line strength S versus magnetic 
field strength at T500 = 0.1 for the 630.2 nm line (Fig. [12^) 
and for the 1564.8 nm line (Fig.[l2b) show that both lines 
are progressively weakened with increasing magnetic field 
strength. The binned averages given by the line indicate 
an almost linear dependence, in the case of the 630.2 nm 
line with some increase of the slope for large field strength. 
The scatter plot of temperature versus field strength (both 
at T500 = 0.1) shown in Fig. [13] indicates that the line 
weakening is largely due to the larger temperature (at 
equ al optical depth) in the m agnetic flux concentrations 
(cf. IChapman fc Sheelev 1968), resulting from the latera. 
radiative heating and reduced density of these structures. 
However, the visible line, which should be more temper- 
ature sensitive owing to its lower value of Xej shows a 
somewhat lower relative weakening than the infrared line. 
This indicates that other effects, such as Zeeman desatu- 
ration as well as vertical velocity and temperature gradi- 
ents, also affect the line strength in a significant way. A 
more detai led consideration on th e basis of response func- 
tions (e.g-. lDel Toro Iniestalliooll Sect. 10.4.2) is required 
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to shed more light on this result, but such an analysis is 
beyond the scope of this paper. 

Fig. [TJ] shows scatter plots of the unsigned Stokes- V 
area, Ay versus the strength of the vertical (line-of-sight) 
component of magnetic field at t 500 = 0.1. It illustrates 
the different magnetic sensitivity of the two lines. The 
Stokes signal of the 630.2 nm line depends linearly on the 
magnetic field strength until the Zeeman splitting is no 
longer small compared to the line width and saturation 
sets in for B z J; 1200 G. The curve for the infrared line 
with its higher Zeeman sensitivity becomes nonlinear at 
lower field strength and bends over around 900 G. This line 
also exhibits a larger scatter of Ay values for intermediate 
field strengths. The decrease of the Stokes- 1^ area of both 
lines for large field strength is due to the temperature- 
induced line weakening in the strong flux concentrations 
(cf. Fig. [T2"]) . The more rapid decrease of Ay with B z 
for the the 1564.8 nm line is consistent with the stronger 
weakening of this line as shown in Fig. 112b . 

Figures [15] and [16] show scatter plots of the Stokes- F 
area asymmetry, 5 Ay, and amplitude asymetry, Say, re- 
spectively, versus the vertical (line-of-sight) magnetic field 
strength. The binned averages of the asymmetries given 
by the red curves show that the asymmetries decrease to- 
wards stronger fields, even with a tendency to become neg- 
ative. The results for the 630.25 nm line are in qualitative 



agreem ent with the observations of Martinez Pille t et al 



(| 19971 see their Figs. 15 and 16) and lSigwarth et alj (]1999l . 
see their Fig. 6). Note, however, the difference in spatial 
resolution between simulation and observations, so that a 
more quantitative comparison is not indicated. 

The asymmetry distributions shown here are in qual- 
itative agreement with the results of Khomenko et al.l 
(20051)) on the basis of a MURaM simulation of a decay- 
ing mixed-polarity field. In both cases (and in agreement 
with observations) it is found that the amplitude asym- 
metry, on average, is a factor 2 to 3 larger than the area 
asymmetry. On the basis of calculat i ons w ith an idealized 
flux tube model. iBellot Rubio et alj (|1997l ) have suggested 
that a downflow in the flux concentration, in addition to 
the external downflow, could lead to asymmetries with 
this property. In fact, the simulated flux concentrations, 
on average, show internal downflows (see Fig. |2T)]) . consis- 
tent with this suggestion. 

The asymmetries for the 1564.8 nm line are systemat- 
ically smaller than those for the visible line. The asym- 
metry is maximized if the wavelength shift due to the 
Zeeman splitting is si milar to the Doppler shift due t o 



the velocity gradient (jGrossmann-Doerth et alj [l988bQ 



This is roughly the case for the visible line while the 
stronger Zeeman sensitivity of the infrared line leads to 
smaller asymmetry compared to that of the visible line. 
Furthermore, the infrared line is weaker and less saturated 
than t he visible line, which additionally reduces the asym- 
metry (|Solankilll989r ). The large scatter of the asymmetry 
values for magnetic field strengths below ~ 200 G is con- 
sistent with the observed as ymmetries of weak Stokes- V 
profiles (e.g.. lSigwarthll2001l ) and indicates that the scat- 



ter in such observations is not entirely due to noise. In 
the simulations, many of the weak V profiles result from 
strongly inclined fields in the granular upflows. Note that, 
according to the definitions of the asymmetries given by 
Eqs. ([3]) and (|4|, complex profiles, e.g. those with 3 or 
more lobes, cannot be distinguished from normal profiles 
in Figs. [15] and [16] Fig. [17] shows that the average am- 
plitude and area asymmetries are well correlated, partic- 
ularly for the 630.2 nm li ne, which is in accorda nce with 
the observational results of Sigwarth et al. ( 19991 see their 
Fig. 9). 

Considering the average F-profile of the 630.25 nm line 
over the whole computational domain of 6 Mm x 6 Mm, 
we find values of 5 Ay = 3.6% for the area asym- 
metry and Say = 17.8% for the amplitude asymme- 
try. This can be compared with similar averages (albeit 
with di fferent averaging area s ) of observed plage spec- 
tra by iMartmez Pillet et all (|l997l . see their Fig. 14), 
who find 5 A v = 3.5 .. .4.5% and day = 9.. .11%. and by 



Sigwarth et al 



19991 . see their Fig. 13), who give values 
of 5 Ay = 3% and Say = 11%. The dominance of posi- 
tive asymmetries (i.e., stronger blue wing of the Stokes- 1^ 
profile) for intermediate field strengths in Figs. [T5l and [TBI 
supports the interpretation that these average asymme- 
tries originate in the peripheral parts of magnetic flux 
concentrations, which expand with height and are sur- 
rounded by strong downflows (such as the example shown 
in Figs. [TCll and [Til see also Table 1). Lines of sight in the 
periphery of a flux concentration cut through the mag- 
netopause and thus exhibit a reversal of gradient of the 
line-of-sight magnetic field (field strength decreasing with 
depth) in comparison to the central parts of the flux con- 
centrations, where the field strength increases with depth 
(see Fig. [T8l) . Together with the increase of the (downflow) 
velocity with depth, this leads to positive asymmetries of 
the Stokes-V^ profiles originating in the peripheral parts 
of the flux concentrations (see the central frame of Fig. [3] 
According to Eq. ©, the sign of the area asymmetry 
should depend on the product of the gradients of magnetic 
field and velocity along the line of sight. Fig. [19] shows 
the dependence of the Stokes-l^ area asymmetry for the 
603.2 nm line on the product (i> z ,r 50 o=o.i — Vz,r soo =o.oi) • 
(I-SzItsoc^o.i _ |Sz|t 500 =o.oi), which corresponds to the av- 
erage gradients over the formation height of the lines, thus 
representing the relevant quantity in Eq. ([8]). It is clear 
that the sign of the average asymmetry is in accordance 
with the prediction from this relation. The scatter results 
from the fact that the quantity considered in Fig. [19] repre- 
sents average gradients over an extended height range and 
thus does not take into account local variations within the 
height range of line formation. The numerous points with 
negative area asymmetry in Fig. [19] in most cases corre- 
spond to very weak Stokes- V profiles formed over granules; 
they give only a minor contribution to the asymmetry of 
the average profile. 

The wavelength shift of the Stokes-F zero cross- 
ing has often been used to determine t he velocity of 
magnetized plasma in observations (e.g., Solankil 19861 : 



S. Shelyag et al.: Stokes diagnostics 



7 



Grossmann-Doerth et al. 119961: iMartmez Pillet et al.. 
1997t ISigwarth et alJll99a) . The 630.2 nm and 1564.8 nm 
lines are formed at somewhat different heights in the solar 
atmosphere, so that they trace the velocity field along 
the line of sight in a different way. The 1564.8 nm line 
is formed deeper, where velocities typically are larger. 
For the simulated magneto-convection, Fig. [20] shows 
scatter plots of the line-of-sight velocity determined from 
the zero-crossing wavelength of the synthetic Stokes-V 
profiles versus the vertical magnetic field at the level 
T500 =0.1. The binned averages given by the line show 
that downflows (negative velocity) dominate for field 
strength above a few hundred Gauss, with the stronger 
flows typically shown by the infrared line. For field 
strength around one kilogauss we find average downflows 
of the order of 1 k rns" 1 with the 630.2 nm line (cf. 
Sigwarth et al. 19991 their Fig. 6) while the deeper orig- 
inating 1564.8 nm shows about 2 kms -1 downflow. The 
weak fields are predominantly associated with granular 
upflows. The strongest redshifts shown by the infrared line 
are at field strengths in the range 500-1000 G, which are 
typical of the expanding field in the periphery of strong 
flux concentrations. For field strengths corresponding to 
the cores of the flux concentrations, the zero-crossing shift 
is reduced to about 1 kms -1 . This suggests that this shift 
could be partly due to the entrainmcnt of material from 
the surrounding granular downflow, an effect probably 
depending also on the spatial resolution of the simulation. 

5. Conclusions 

We have analyzed the spectro-polarimetric signatures of 
3D MHD models of solar magneto-convection, based on 
the synthesis of Zeeman-sensitive Fel lines in the visi- 
ble and in the infrared. The comparison between various 
properties of the synthetic Stokes profiles (such as line 
strength, line width, Doppler shift, Stokes-V signal, area 
and amplitude asymmetries) with the physical properties 
of the model (temperature, magnetic field, and velocities) 
provides insight into the mechanisms underlying various 
characteristics of the observed spectra. In particular, the 
generation of strongly asymmetric Stokes-V profiles in the 
periphery of magnetic flux concentrations largely confirms 
earlier simplified models based upon lines of sight cross- 
ing the magnetopause of flux concentrations fanning out 
with height. The line weakening in strong magnetic fields 
is found to be related to larger temperature (and smaller 
temperature gradient) in the simulated magnetic flux con- 
centrations, in agreement with earlier interpretations of 
observational results. 

Acknowledgements. The authors are grateful to Sanja 
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Fig. 1. Vertical magnetic field component (left panels) and vertical velocity (right panels) on the surface T500 = 1 
(continuum optical depth unity at 500 nm wavelength) for snapshots from simulations with (B z ) = 200 G (plage, 
upper row) and (B z ) — 10 G (decay phase, bottom row). Positive values of the velocity (shown in blue) correspond to 
upflows; downflows are are indicated in red. Note the different scales for positive and negative B z in the left panel (all 
negative-polarity field is rather weak) . 
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Fig. 2. Comparison between the spatially averaged syn- 
thetic Fei line profiles for the simulation with (B z ) — 10 G 
(diamond symbols) and observed quiet-Sun profiles (soli d 
curves) from the spectral atlas of iDelbouille et alJ (|1973l ). 
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Fig. 3. Properties of synthetic Stokes-/ and Stokes- ^profiles for the snapshot from the (B z ) — 200 G simulation shown 
in Fig. [TJ Given are maps of the (normalized) continuum intensity at 630.2 nm (top left panel), line strength, S (top 
middle panel for Fei 630.2 nm and top right panel for Fei 1564.8 nm), unsigned Stokes- area, Ay (middle left panel 
for 630.2 nm and bottom left panel for 1564.8 nm), Stokes- area asymmetry (central panel for 630.2 nm and bottom 
middle panel for 1564.8 nm), and Stokes- y amplitude asymmetry (middle right panel for 630.2 nm and bottom right 
panel for 1564.8 nm). In the continuum intensity image (top left panel), the slit position for Figs. [5] and is indicated 
by the vertical line and the cut through the magnetic flux sheet shown in Fig. [TO] is indicated by the inclined line in 
the upper left quadrant. 
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Fig. 4. Same as Fig. [3] for the snapshot from the simulation with (B z ) — 10 G. The scales are the same as in the 
previous figure, except for the Ay images. 
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Fig. 5. Synthetic slit spectra, calculated for the 630.2 nm line. Top left: Stokes-/, top right: Stokes- V/I, bottom left: 
Stokes-Q/J, bottom right: Stokes-C//J. The position of the slit is indicated by the vertical line in the top left panel of 
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Fig. 6. Same as Fig. [5] for the 1564.8 nm line. 
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Fig. 7. Vertical velocity (top left panel) and field strength (top right panel) at the surface T500 = 0.1 for the (B z ) — 
200 G simulation. The middle row shows maps of the vertical velocity as determined from the Doppler shift of the 
630.2 nm line (middle left panel) and the 1564.8 nm line (middle right panel). The bottom row gives maps of the line 
width of the 630.2 nm line (bottom left panel) and the 1564.8 nm line (bottom right panel). 
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Fig. 8. Same as Fig. [7] for the (B z ) = 10 G simulation. 
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Fig. 9. Scatter plots of velocity determined from the Doppler shift of the center of gravity of the spectral line vert 
the vertical velocity at T500 = 0.1 from the simulation (200 G). The line shows the binned average. 
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Fig. 10. Temperature (upper left panel), gas pressure (upper right panel), magnetic field strength (lower left panel), 
and vertical velocity (lower right panel) for the vertical cut through the simulation box marked by the white line on 
the upper left panel of Fig. [3l The green curve indicates the level T500 = 1- Negative values of the velocity (shown in 
red) correspond to downflows. The lines of sight corresponding to the Stokes profiles and gradients shown in Fig. [TT] 
are indicated at the top of each panel. 
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Fig. 11. Stokes-/ and Stokes- V profiles for the 1564.8 nm line (third and fourth column) and for the 630.2 nm line 
(fifth and sixth column). The corresponding magnetic field and velocity profiles as a function of logrsoo are given in 
the first two columns. The level T500 = 1 is marked by the dash-dotted line; zero vertical velocity is indicated by the 
dashed line. Each row corresponds to one of the vertical lines of sight indicated in Fig. [TU1 ranging from the central 
part of the flux concentration (top row) to its outer periphery (bottom row). 
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Fig. 12. Scatter plots of line strength versus magnetic field strength at T500 = O.f for the snapshot from the (B z ) = 
200 G simulation, a: Fei 630.2 nm, b: Fei 1564.8 nm. The lines show the binned averages. Note the different vertical 
scales on the two panels. 
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Fig. 13. Scatter plots of temperature versus magnetic field strength (both at T500 = 0.1). The line shows the binned 
average. 
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Fig. 14. Scatter plots of the unsigned Stokes- V area versus the vertical (line-of-sight) component of the magnetic field 
at T500 = 0.1 a: Fei 630.2 nm, b: Fel 1564.8 nm. The lines show binned averages. 




Fig. 15. Scatter plots of the Stokes-V area asymmetry versus the vertical (line-of-sight) component of the magnetic 
field at T500 = 0.1 a: Fei 630.2 nm, b: Fel 1564.8 nm. The lines show binned averages. 




Fig. 16. Same as Fig. [15] for the Stokes- V amplitude asymmetry. 
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Fig. 18. Spatial distribution of 

(B 2iT500=0 .i ~ B z . Two=om ) 1 representing the average 
line-of-sight gradient of the vertical magnetic field. 
Magnetic field increasing with optical depth is indicated 
by blue color, magnetic field decreasing with optical 
depth is shown in red. 
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Fig. 19. Scatter plot of the Stokes- area asymmetry for 
the 630.2 nm line versus a quantity representing the gradi- 
ents of the line-of-sight magnetic field and velocity within 
the height range of line formation. The sign of the bin- 
averaged 5 Ay (indicated by the line) is in agreement with 
the relation given in Eq. (JSJ). 
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Fig. 20. Scatter plots of the line-of-sight velocity calculated from the zero-crossing-wavelength of the Stokes- V profile 
versus the vertical component of the magnetic field at T500 = 0.1. a: Fei 630.2 nm, b: Fei 1564.8 nm. The lines show 
binned averages. On average, strong magnetic flux concentrations are associated with downflows (negative velocity), 
while weak fields B < 200 G show some preference for (granular) upflows. 



